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Abstract 


We  present  a  new  approach  for  inference  from  accelerated  life  tests. 
Our  approach  is  based  on  a  dynamic  general  linear  model  setup  which  arises 
naturally  from  the  accelerated  life  testing  problem  and  uses  linear  Bayesian 
methods  for  inference-  The  advantage  of  the  procedure  is  that  it  does  not 
require  large  number  of  items  to  be  tested  and  that  it  can  deal  with  both 
censored  and  uncensored  data.  Furthermore,  the  approach  produces  closed 
form  inference  results.  We  illustrate  the  use  of  our  approach  with  some 

actual  accelerated  life  test  data.  - _____ - 
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1.  INTRODUCTION  AND  OVERVIEW 

In  reliability  studies  it  is  a  common  practice  to  subject  items  to  an 
environment  which  is  more  severe  than  the  normal  operating  enviroment  so 
that  failures  can  be  induced  in  a  short  amount  of  test  time.  A  more  severe 
environment  can  be  created  by  increasing  one  or  more  of  the  stress  levels, 
which  constitute  the  environment,  to  values  which  are  greater  than  their 
usual  levels.  Such  tests  are  called  accelerated  life  tests.  The  main  problem 
with  inference  from  accelerated  life  tests  is  that  uncertainty  statements 
about  the  failure  behavior  of  the  items  at  usual  stress  conditions  have  to  be 
made  using  life  length  data  from  the  more  severe  stress  conditions.  Most  of 
the  current  literature  on  the  accelerated  life  testing  problem  is  based  on  the 
sample  theoretic  paradigm  [see  Meeker  and  Hahn  (1985)  for  an  up-to-date 
review]  ;  exceptions  to  these  are  Meinhold  and  Singpurwalla  (1984)  and 
Blackwell  and  Singpurwalla  (1986)  which  involve  a  Kalman  Filter  formulation 

of  the  accelerated  life  testing  problem. 

In  this  paper  we  present  a  new  procedure  for  inference  from 
accelerated  life  tests.  Our  procedure  is  Bayesian  and  is  based  on  a  dynamic 
general  linear  mode!  (DGLM)  setup  [see  West,  Harnson  and  Migon  (1985), 
henceforth  WHM]  which  arises  naturally  from  the  accelerated  life  testing 
scenario.  Our  procedure  uses  the  linear  Bayesian  approach  of  WHM  and 
produces  closed  form  solutions  for  inference.  The  main  advantage  of  the 
procedure  is  that  it  does  not  require  large  number  of  items  to  be  tested  at 
each  stress  level  and  that  it  can  deal  with  both  censored  and  uncensored 
data.  Besides,  the  recursive  nature  of  the  results  produced  by  the 
procedure  facilitates  its  use  on  Personal  Computers. 
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In  our  setup  we  assume  that  lifetimes  obtained  under  all  stress  levels 
have  an  exponential  distribution  and  make  use  of  a  particular  time 
transformation  function.  The  extension  of  our  approach  to  other  failure 
distributions  such  as  the  Weibull  distribution  is  being  investigated.  The 
extension  of  the  procedure  to  other  time  transformation  functions  is  straight 

forward. 

A  synopsis  of  our  paper  is  as  follows: 

In  Section  2  we  present  the  notation  and  preliminaries  for  the 
accelerated  life  testing  problem.  In  Section  3  we  introduce  the  power  law  as 
a  time  transformation  function  and  present  the  OGLM  setup  for  the  problem. 
We  also  discuss  the  adoption  of  our  procedure  to  other  time  transformation 
functions.  In  Section  4  we  describe  the  inference  procedure  for  the  DGLM 
setup.  Finally  in  Section  5,  we  illustrate  the  use  of  our  approach  by 
applying  it  to  some  actual  accelerated  life  testing  data  published  by 
Nelson  (1972). 

2.  NOTATION  AND  PRELIMINARIES 

Assume  that  testing  is  done  in  stages  using  k  accelerated  stress 
levels  of  decreasing  intensity  which  are  specified  in  advance  and  may  or  may 
not  be  equally  spaced.  During  the  i-th  stage  of  testing,  items  are  tested 
under  an  accelerated  stress  level  denoted  by  S^.  Using  the  notation  S<  <  S^ 
to  denote  that  S^  is  a  more  severe  testing  environment  than  S,,  we  note  that 

Si  >  S2  >  -•  >  >  Sfc+i  ,  (2.1 

where  S»_+i  denotes  the  use  stress  at  which  no  testing  is  done. 
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At  each  stress  level  S„  n,  items  are  tested  for  a  predetermined  and 
fixed  time  length  t,.  Let  Y„  denote  the  time  to  failure  of  the  j-th  item 
tested  under  stress  level  S„  and  denote  the  realization  of 
j-Qj.  j_  2,  ,  r,  <  n^  where  r^  denotes  the  number  of  failures  observed 

during  the  time  interval  (0,  T^l. 

We  assume  that  the  failure  rate  function  for  items  tested  under 
stress  S,  is  constant  and  denoted  by  Xj.  Thus,  given  X,  the  failure 
distribution  under  stress  S,  is  described  by  the  exponential  density 

p(yslXt)  =  x,e“  . 


Furthermore,  given  X^,  failure  times  for  items  tested  under  stress  Sj  are 
judged  to  be  independent. 

Given  the  rj’s  and  yj/s,  our  goal  is  to  make  inferences  about  the 
failure  behavior  of  an  item  operating  at  use  stress  (normal)  conditions  S^+i. 

In  so  doing,  it  is  most  common  to  assume  a  functional  relationship  between 
the  failure  rate  and  the  applied  stress  level.  Such  relationship  is  known  as 
an  acceleration  or  time  transformation  function.  Commonly  used  models  for 
describing  such  relationship  are  the  Arrhenius  Law.  the  Eyring  Law,  and  the 
Power  /aw;  see  for  example  Mann,  Schafer  and  Singpurwalla  (1974),  p.  421. 

In  what  follows,  we  will  focus  attention  on  the  popular  Power  Law  and  present 
the  DGLM  setup  for  the  accelerated  life  testing  problem.  In  general,  the  use 
of  any  of  these  laws  should  be  based  on  the  physics  of  failure  for  the 
problem  at  hand. 
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3.  THE  DGLM  SETUP  FOR  ACCELERATED  LIFE  TESTING 


Under  the  Power  Law,  we  write 


—  » 


(3.1) 


where  a  and  /?  are  unknown  coefficients  which  describe  the  stress  effect  on 
failure  rate.  The  subscripts  associated  with  a  and  0  imply  the  fact  that  the 
time  transformation  function  might  be  changing  from  one  stress  level  to 
another.  As  noted  by  Blackwell  and  Singpurwalla  (1986)  this  is  likely  to 
happen  due  to  the  changes  in  the  basic  failure  mechanism  with  changes  in  the 
stress  level. 

We  can  linearize  (3.1)  by  taking  natural  logarithms  on  both  sides  and 

write 


r),  =  loga,  +  0t  logS<  , 


(3.2) 


where  T)i  =  logX<  .  We  define  (loga^  0,)  ,  F.t'=  d  logSj)  and  write 

(3.2)  as 


We  note  that  (3.3)  provides  us  with  a  guide  relationship  in  the  sense  of  WHM 
and  that  is  the  underlying  state  vector. 

^gxt  we  describe  how  the  time  transformation  function  is  changing 
from  one  stress  level  to  another  by  specifying  the  system  (evolution) 
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equation  of  the  model  as 


,  (3.4) 

where  Wj  is  a  random  innovation  term.  We  note  that  in  (3.4)  the  state  vector 
6j’s  are  assumed  to  be  constant  from  one  stress  level  to  another,  except  for 
the  random  changes  brought  about  by  the  innovations  The  distribution 
of  w,  is  only  partially  specified  through  its  first  and  second  order  moments. 
We  use  the  notation. 


w«  ~  I  0  .  ]  ,  (3.5) 

to  denote  the  fact  that  jwj  has  mean  vector  0  and  variance-covariance  matrix 
W{.  Typically,  jwj  is  independent  of 

Let  D,  denote  all  the  relevant  information  available  at  stage  i  after 
observing  r,  and  failure  times  {yt^}^  ^  at  stress  level  S,;  Do  represents  all 
relevant  information  available  at  stage  0  prior  to  any  testing.  At  stage  (i-1), 
we  assume  that  the  distribution  of  the  state  vector  is  partially 
described  by  the  first  and  second  order  moments  as 

( ^{—1 1  (^1—1  )  I  J5,t— 1  »  £»— 1  ]  •  (3.6) 

Using  the  system  equation  and  (3.6)  we  can  write 

(  )  ~  (  m,_i  ,  Rt  ]  ,  (3.7) 
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where  Rj  =  C,_,  +  W,  .  Equation  (3.7)  provides  us  with  a  partial 
description  of  the  prior  distribution  of  the  state  vector  before  testing  at 
stress  level  S,.  It  is  important  to  note  that  the  innovation  term  w,  provides 
an  increase  in  uncertainty  (represented  by  the  addition  of  Wf)  over  the 
stress  levels  as  changes  to  This  loss  of  information  about 
motivates  the  discount  concept  (see  for  example,  Smith  (1979)  and  Ameen  and 
Harrison(1985)]  used  by  WHM  as  a  guide  for  the  choice  of  W,.  The 
underlying  idea  is  that  the  increase  in  uncertainty  over  the  next  stress  level 
should  be  relative  to  that  available  at  the  present  stress  level,  measured  by 
C,_i.  The  details  associated  with  the  choice  of  the  discount  factor  are 
described  by  WHM.  In  our  case,  the  discount  factor  can  be  chosen  as 
function  of  (St/S,_i)  to  take  into  account  the  relative  magnitudes  of  the 

stresses. 

From  the  guide  relationship  (3.3)  we  can  write 


fj  =  ElJ7jD,_i]  =  F/m4_i  , 

Qi  =  Var(^{|Dj_i]  =  ' 


(3.8) 


We  note  that  (3.8)  provides  us  with  only  the  first  two  moments  of  the  prior 
distribution  of  Vi  at  stage  (i-1).  We  can  specify  a  full  distributional  form  for 
the  prior  of  Tfi  and  use  the  relationships  given  by  (3.8)  to  determine  the 
parameters  of  the  distribution.  In  order  to  obtain  closed  form  solutions  in 
the  Bayesian  analysis,  we  will  specify  a  conjugate  form  for  the  prior  of  Vi 
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which  is  a  loggamtna  density  of  the  form 


p(77,  ID,_,)  oc  exp{  a,?7,-  b,e^‘  },  (3 

where  a^  and  bj  are  the  prior  parameters.  We  denote  the  density  in  (3.9)  by 
JtQ(a,  ,bi).  From  (3.9)  we  may  obtain 

E[r7, 1  D,_i]  =  ♦(a,)  -  log(b*), 

and 

Var(r7jD,_,l  =  ♦'(aa 


where  ♦(  -)  and  ♦'(  •)  are  the  dIgamma  and  trlgamma  functions  respectively. 
We  will  specify  the  prior  parameters  a<  and  b,  such  that  the  first  two 
moments  of  agree  with  (3.8),  that  is. 


♦(at)  —  log(bt)  =  f„ 

♦'(at)  =  q,. 


(3.10) 


In  solving  (3.10)  for  aj  and  bt  some  approximations  to  ♦(  • )  and  ♦'(  )  can  be 
used  [see  for  example,  Cox  and  Lewis  (1966),  Ch.  2].  We  note  that  specifying 
the  prior  distribution  of  Vi  as  JL<3(at,bt)  implies  that  the  prior  of  Xt  at  stage 
(i-1)  is  a  gamma  density  with  shape  parameter  at  and  scale  parameter  bj.  We 
will  denote  this  density  as  Q(at,bt). 
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To  summarize,  the  DGLM  setup  for  the  accelerated  life  testing 


problem  is  as  follows: 


p(y£  IM 


Vi 
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(rjJDi.i)  ~  Vt  =  logCM,  (3.11) 

(§,1  ID<— 1'  f  5.*— 1  ’  .5.<  ^ 


where  a(  and  bj  are  chosen  according  to  the  guide  relationship  TJi  —  jF/  so 

that  E[r7i|D,_i]  =  fi  and  Var[r?,  jD<_i]  =  Qj. 

We  note  that  the  DGLM  setup  for  the  accelerated  life  testing  problem 
can  be  easily  modified  for  other  time  transformation  functions.  For  example, 
under  the  Arrhenius  Law 


=  exp{ai— ^i/SJ  ,  (3.12) 

which  implies  that  the  guidance  relationship  =  log(K^)  is  given  by 
equation  (3.3)  with£'=  (a,  00  and  F/=  (1  -1 /S^).  Now  the  DGLM  setup 
of  (3.11)  can  be  applied  to  the  problem.  Extension  to  the  Eyring  and  other 
laws  follows  along  the  same  lines. 

4.  INFERENCE  RESULTS  FOR  DGLM  SETUP 

Given  D„  the  available  information  after  the  i-th  stage  of  testing,  we 
need  to  update  our  inferences  about  T]^  and  The  posterior  distribution  of 
TJi  given  Dj  is  obtained  via  the  standard  use  of  Bayes'  theorem 

p(r7jDt)  a  m;D<)  p(r7tlDj_i), 
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where  Z.f?7,;Di)  is  the  likelihood  function  for  given  by 


I(r?t;D,)  =  exp{r7,r,  -  T,e^‘}. 


(4.1) 


and  T,  =  2  Vu  +  stress  level  S,.  It 

j-i 

can  be  shown  that 


(riJD.)  ~  X.g(at+r,,b,+T,), 


(4.2) 


implying  that  (XJD()  ~  Q(aj+rt,b, +T,).  The  posterior  mean  and  variance  of 
T]^  are  given  by 


g,  =  E(r?,|D,]  =  ♦(a,+r.)  -  log(b,+T,), 


(4.3) 


p,  =  Varlr/,|DJ  =  ♦'(a<+r,). 


The  next  step  is  to  obtain  the  posterior  mean  and  variance  of  the 
state  vector  0{.  We  note  that  the  full  form  of  the  posterior  distribution  for 
(6,  ID,)  is  not  available  since  the  prior  of  (ejDt-i)  is  only  partially  specified. 
Recognizing  the  fact  that,  given  TJt,  the  observation  model  (2.2)  does  not 
depend  on  WHM  developed  a  method  for  updating  the  first  two  moments 
of  (6,  ID,)  using  the  linear  Bayesian  approach  of  Hartigan  (1969).  The  details 
of  the  method  is  omitted  here  and  only  the  main  result  is  presented. 


9 


The  method  described  by  WHM  leads  to  the  updating  of  the  state 


vector  as 


where 


(0,  ID,)  ~  [m,  ,  C,|, 


(4.4) 


m,  =  m,_,  +  fi)/qi» 


q,'  ’ 


(4.5) 


with  s,  =  iR,  F,.  We  note  that  the  updating  equations  for  in  the  DGLM 
setup  have  the  same  form  as  the  well-known  Kalman  filter  recursions  [see  for 
example  Meinhold  and  Singpurwalla  (1983)]. 

After  performing  the  test  at  the  last  stress  level  S^,  our  aim  is  to 
make  inference  about  the  life  length  of  an  item  at  use  stress  conditions  S;t+i. 
In  other  words,  given  we  wish  to  make  uncertainty  statements  about 
and  Xfc+i  (or  equivalently  about  Vk+0’  Given  Dit  we  have 

(r/JD^)  ~  Xg(ak+rk,bk+T;t), 

and 

(  f  Sk  »  £k 


We  describe  our  uncertainty  about  the  state  vector  using  the  system 

equation;  that  is. 


IDjt)  ~  (  jEfc  9  ^fc+1  ]»  ^ 

where  R^+j  =  C*.  +  From  the  guide  relationship  of  (3.3),  we  obtain 
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=  F'k+1  and  q-t+i  ==  ■  Using  the  DGLM  setup  of 

(3.11),  we  describe  uncertainty  about  Vk+i  density 

(^k  +  llUjt)  —  ^Ql®k  +  1  •  ^k+ll* 


where  a^+i  and  b^+i  are  chosen  such  that  ♦(at+i)  —  log(b^+i)  =  f^+i 
and  ♦'(afc+i)  =  q^+i-  Thus,  inference  about  the  use  stress  failure  rate  is 
made  by  using  the  probability  density  (^it+iIDj.)  Q(®t+i»  bk+il- 

Inference  about  the  life  length  of  an  item  operating  at  use  stress 
conditions  is  obtained  using  the  predictive  density 

oo 

p(yik+ilU]k)  =  I  p(yk+i  l^k+i)  p(^k+ilF>k)  dXfc+i 

0  (4.8) 

_ ®k+l(^k+l) 

(l>k+l+  yk+l) 


which  is  a  Pareto  density  with  reliability  function 


R(yk+ilUfc)  — 


(b^+i) 


*k+l 


(bit+i4-  yk+i) 


»k+i 


(4.9) 


A  nice  feature  of  the  DGLM  setup  is  that,  at  any  stress  level  Sj,  we 
can  make  inference  about  the  failure  behavior  at  any  stress  <  S<  (  that  is, 
i  <  j  )  and  at  the  use  stress.  Assume  that  we  have  tested  items  at  stress 
level  St  and  we  have  (€-1)  stages  of  testing  ahead,  that  is,  £  =  k— i  +  1.  We 
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define 


and 

C,(£)  =  Var(£^^^|D.], 

where  6,(0)  =  m,  and  C,(0)  =  C,.  Using  the  system  equation  of  the  DGLM 
setup  we  obtain 


mi(£)  =  m,  ,  for  £  ^  0  , 

C,(£)  =  £,(£-!)  +  W  for  £  >  0. 

Similarly,  using  the  guide  relationship  we  write 


(4.10) 


f,(£)  =  E{  77.,^^  ID,  ]  =  *.(«). 

q,(£)  =  Var[  77^^^|D.  ]  =  £ ',+£ £,(£)  F^^£. 


(4.  11) 


Thus,  we  can  use  (4.10)  and  (4.11)  in  the  DGLM  setup  and  make 
inference  about  the  failure  rate  at  stress  levels  Sj  <  S,  given  D,. 


5.  EXAMPLE 

The  method  was  applied  to  some  accelerated  life  test  data  taken  from 
Nelson  (1972).  This  data,  given  in  Table  5.1,  represents  the  times  to 
breakdown  of  an  insulating  fluid  with  subjected  to  various  voltage  levels. 
The  insulating  fluid  is  tested  at  accelerated  stress  levels  26,  28,  30,  32,  34, 
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36,  and  38  Kv,  and  inference  is  to  be  made  concerning  the  breakdown  times 
for  the  the  insulating  fluid  at  20  Kv  (use  stress).  In  our  approach  we  assume 
that  testing  was  done  in  a  sequential  manner  from  highest  to  lowest  stress 
level. 

TABLE  5.1 

Times  to  Breakdown  of  an  Insulating  Fluid  (in  Minutes) 

Under  Various  Values  of  the  Stress 


38Kv 

36Kv 

34Kv 

32Kv 

30Kv 

28Kv 

26Kv 

.09 

.35 

.19 

.27 

7.74 

6835 

5.79 

.39 

.59 

.78 

.40 

17.05 

108.29 

157932 

.47 

.96 

.% 

.69 

20.46 

110.59 

2323.70 

.73 

.99 

1.31 

.79 

21.02 

426.07 

.lA, 

1.69 

2.78 

2.75 

22.66 

1067.60 

1.13 

1.97 

3.16 

3.91 

43,40 

1.40 

2.07 

4.15 

9.88 

47.30 

2.38 

2.59 

4.67 

13.95 

139.07 

2.71 

4.85 

15.93 

141.12 

2.90 

6.50 

27.80 

175.88 

3.67 

735 

5334 

194.90 

3,99 

8.01 

8235 

5.35 

8.27 

89.29 

13.77 

12.06 

100.58 

25.50 

31.75 

215.10 

3232 

33.91 

36.71 

7239 

Due  to  our  unfamiliarity  with  the  problem  at  band,  the  selection  of  the  prior 
parameter  values  is  made  for  illustrative  purposes  only.  We  use  the  Power  Law  as 
our  time  transformation  function  and  specify  our  prior  parameter  values  as 

0.5  0.0 

0.0  0.5 


=  |^-15.0  10.0  j. 


Co 
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and 

0.01  0.0 
0.0  0.01 

for  all  stages  i. 

Plots  of  the  posterior  means  of  Log  a,  and  (which  are  the  components  of  m^) 
for  each  stage  of  testing  are  presented  in  Figures  5.1  and  5.2  respectively.  Note 
that  there  is  an  overall  downward  trend  in  both  graphs  and  that  after  a  strong 
decline  from  the  prior  mean  values,  a  more  stable  behavior  is  exhibited.  Thus 


nsTiNc  sma 

Fig.  5.1  Posterior  Means  of  log  at 


implying  that  our  prior  belief  about  the  effect  of  stress  on  failure  behavior  of  the 
units  was  stronger  than  that  exhibited  by  the  data. 

In  Table  5.2  we  present  the  E(K^|  DJ  ,  j  =  i,  i+l»  •••  »  L,  k  +  l  for  each  testing 
stage  i.  Columns  of  this  table  exhibit  our  assessment  of  the  failure  behavior  of  the 
system  (insulating  fluid)  under  the  different  levels  of  stress  based  on  the  total 
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Fig.  5.2  Posterior  Means  of 


amount  of  information  available  after  each  stage  of  testing.  In  Figure  5.3  we  plot  the 
last  row  of  this  table  to  illustrate  the  evolution  of  the  failure  rate  predictions  for 
use  stress.  In  addition,  in  Figure  5.4  we  illustrate  the  predictive  reliability  functions 
for  use  stress  after  testing  stages  1,  3,  5,and  7. 


TABLE  5.2 

Failure  Rate  Predictions  at  Different  Testing  Stages 


PREDICTION 
FOR  STAGE 

1 

2 

TESTING 

3 

STAGE 

4 

5 

6 

7 

1 

1.1102 

2 

0.8155 

0.2546 

3 

0.6288 

0.1944 

0.0810 

4 

0.4773 

0.1503 

0.0622 

0.0293 

5 

0.3559 

0.1142 

0.0480 

0.0234 

0.0154 

6 

0.2601 

0.0852 

0.0363 

0.0171 

0.0115 

0.0047 

7 

0.1857 

0.0622 

0.0270 

0.0129 

0.0087 

0.0034 

0.0015 

8  (use 

0.0563 

0.0204 

0.0094 

0.0047 

0.0033 

0.0014 

0.0006 

15 


6 


ACKNOWLEDGEMENTS 


This  rssearch  was  supported  by  Contract  N00014-85-K0202,  Project  NR  042-372, 
Office  of  Naval  Research  and  by  Grant  DAAL  03-87-K-0056,  U.  S.  Army  Research 
Office, 


REFERENCES 

[1]  J.  R.  M.  Ameen,  P.  J.  Harrison,  “Normal  discount  Bayesian  Models”, 

in  Bayesian  Statistics  2,  eds.  J.  M.  Bernardo,  M.  H.  DeGroot,  D.  V.  Lindley, 

A.  F.  M,  Smith,  North— Holland,  1985,  pp.271-298  . 

[2]  L.  M.  Blackwell,  N.  D  Singpurwalla,  “Inference  from  accelerated  life  tests  using 
filtering  in  colored  noise”.  Technical  Report  GWU/lRRA/TR-86/11,  Institute  for 
Reliability  and  Risk  Analysis,  George  Washington  University,  1986. 

[3]  D.  R.  Cox,  P.  A.  W.  Lewis,  Statistical  Analysis  of  Series  of  Events, 

Methuen,  1%6. 

[41  J.  A.  Hartigan,  “Linear  Bayesian  Methods”,  Journal  of  the  Royal  Statistical 
Society,  Series  B,  vol  31,  1%9,  pp  446-454. 

[5]  N,  R.  Mann,  R.  E.  Schafer,  N.  D.  Singpurwalla,  Methods  for  Statistical  Analysis 
of  Reliability  and  Life  Data,  John  Wiley  and  Sons,  1974. 

[6]  W.  Q.  Meeker,  G.  J.  Hahn,  “How  to  plan  an  accelerated  life  test  —  some  practical 
guidelines”,  American  Society  for  Quality  Control,  vol.  10,  1985, 

[7]  R.  E.  Meinhold,  N.  D.  Singpurwalla,  “Understanding  the  Kalman  filter”, 
applications”,  American  Statistician,  vol  37,  1983,  pp.  123-127. 


17 


[8J  R.  E.  Meinhold,  N.  D.  Singpurwalla,  “A  Kalman  filter  approach  to  accelerated  life 
testing  —  a  preliminary  development”,  in  Stochastic  Failure  Models, 
Replacement  and  Maintenance  Policies,  and  Accelerated  Life  Testing, 
ed.  M.  Abdel-Haneed,  E.  Cinlar,  J.  Quinn,  Academic  Press,  1984,  pp.  169-175. 

[9]  W.  B.  Nelson,  “Graphical  analysis  of  accelerated  life  test  data  with  the  inverse 
power  law  model”,  IEEE  Trans.  Reliability,  R-21,  1972,  pp.2-11. 

110]  J.  Q.  Smith,  “A  generalization  of  the  Bayesian  steady  forecasting  model”. 

Journal  of  the  Royal  Statistical  Society,  Series  B,  vol  41,  1979,  pp  3JS-3S7. 

(11]  M.  West,  P.  J.  Harrison,  H.  S.  Migon,  “Dynamic  generalized  linear  models  and 
Bayesian  forecasting”.  Journal  of  the  American  Statistical  Association, 
vol  80,  1985,  pp  13-91. 


18 


